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ABSTRACT 



We show that the starburst-driven outflow in the peculiar galaxy NGC 2782 forms a 
well-defined coUimated bubble which has an extent of ~ 1 kpc and a closed shell at its edge, 
as seen in Ha, [O III], and 5 GHz radio continuum. The shell coincides with the maximum in 
intensity and linewidths of [O III] lines in a blueshifted emission nebula which was previously 
detected via optical spectroscopy by Boer ct al. (1992). Such a remarkable outflow morphology 
has not been observed to date in any other starburst galaxy of comparable luminosity. The 
radio continuum map reveals a second bubble of similar size on the opposite side of the nucleus, 
forming a striking double-bubble outflow morphology. We argue from the morphology and short 
timescale (~ 4 x 10^ years) of the outflow that it is dynamically younger than freely-expanding 
outflows seen in other galaxies which harbor circumnuclear starbursts of comparable luminosity, 
e.g., M82. We suggest that the outflow in NGC 2782 is in the early stage where thermal 
instabilities have not yet completely ruptured the outflow bubble. We present evidence that 
the outflow is driving warm and hot ionised gas, and possibly cold molecular gas, out of the 
central kpc of the galaxy. We estimate the contribution of the hot, warm, and cold phases of 
the ISM to the energetics of the outflow. This study is based on our optical BVR, Ha, and 
[O III] observations from the WIYN telescope and OVRO CO interferometric data, along with 
available 5 GHz radio continuum and ROSAT X-ray maps. 

Subject headings: galaxies: starburst galaxies: ISM galaxies: interactions — 
ISM: jets and outflows — galaxies: evolution — galaxies: structure 



- 3 - 



1. Introduction 

St arburst- driven winds can influence the formation and evolution of galaxies. Such winds 
can drive mctal-cnriched material out of the galactic plane, thus affecting the mass-mctallicity 
relation between galaxies, the metallicity-radius relation within galaxies and the enrichment of 
intracluster medium (e.g, Dekel & Silk 1986; Heckman, Lehnert, & Armus 1993, and references 
therein). In addition, the properties of young stellar components built by central starbursts 
depend on the rate at which starburst-driven winds blow the ISM out of the galaxy. It is therefore 
important to study how starburst-driven winds are triggered and evolve. A starburst-driven 
outflow can be associated with a large spectrum of emission e.g., radio continuum (RC), optical, 
UV, and X-rays. The relative intensity and distribution of the emission at different wavelengths 
give critical clues about the evolutionary phase of the outflow and physical conditions in the 
ISM. Nuclear outflows have been observed in a large number of starburst and Seyfert galaxies 
(e.g., M82, NGC 253, NGC 3628, NGC 3079, Arp 220, NGC 6240, NGC 4666, NGC 4945, 
NGC 1068, NGC 1320, NGC 4235) as discussed by many authors (see Heckman et al. 1993; 
Baum et al. 1993; Lehnert & Heckman 1995; Colbert et al. 1996; Dahlem et al. 1997, & 
references therein). In this paper, we describe the starburst-driven outflow in the peculiar galaxy 
NGC 2782; the outflow has a remarkable morphology that has not been observed to date in any 
other starburst galaxy of comparable luminosity. This study is based on our optical BVR, Ha, 
and [O III] images from the WIYN (Wisconsin Indiana Yale NOAO) telescope and 2" resolution 
OVRO (Owens Valley Radio Observatory) CO (J=l->0) data, along with available 5 GHz RC 
observations (Saikia et al. 1994), and an archival ROSAT X-ray map. In this paper, we focus on 
the starburst-driven outflow and describe the optical observations in detail. In paper II (Jogee, 
Kenney & Smith 1997), we describe the observations and full analysis of the CO data, and 
present evidence for a nuclear bar fuelling molecular gas into the central starburst. 

2. Previous work on the inner few kpc of NGC 2782 

NGC 2782 is a nearby (D=34 Mpc for Hq of 75 km Mpc^^), peculiar galaxy with a 
pair of HI and optical tails (Sandage & Bedke 1994; Smith 1991), and three ripples within the 
disk (Arp 1966; Smith 1994). The HI and optical properties can be modelled with a recent 
interaction or merger of two disk galaxies of unequal mass (Smith 1994). Previous Ha maps 
(Hodge & Kennicutt 1983; Smith 1994; Evans et al. 1996) showed an arc of emission along the 
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inner ripple at r ~ 25 " (4.2 kpc), and an unresolved bright region of star formation in the inner 
10" (1.7 kpc) radius. This region harbors one of the most luminous circumnuclear starbursts 
among nearby (D < 40 Mpc) spirals (Devereux 1989), with a FIR luminosity (2x IO^^Lq; Smith 
1991) comparable to M82. The circumnuclear region has optical emission spectra indicative 
of HII regions (Sakka et al. 1973; Balzano 1983; Kinney et al. 1984), as well as an additional 
component of highly excited, highly ionised gas (Kennicutt et al. 1989). Boer, Shultz, and Keel 
(1992) used spatially resolved spectra to show that the high excitation gas lies in a bhicshifted 
outflowing emission nebula. The 5 GHz RC map at 1" resolution presented by Saikia et al. 
(1994) revealed an intriguing set of peaks and bubbles, for which no explanation has yet been 
proposed. The detailed distribution and kinematics of molecular gas in the central region was 
previously unknown as NGC 2782 had been mapped only with single dish observations (e.g.. 
Young et al. 1995) and low resolution (6") interferometry with the Nobeyama Millimeter Array 
(Ishizuki 1994). The latter study shows an elongated feature with two barely resolved CO 
peaks. The higher (2") resolution CO data we present in paper II (Jogee et al. 1997) resolves 
the elongated feature into a clumpy, double-lobed bar-like feature of radius ^ 7.5" (1.3 kpc). 

3. The remarkable outflow morphology revealed by new observations 

Two 5 min. exposures of NGC 2782 were taken in Ha+ [N II] , [O III] , off-line continuum 
and Harris BVR filters on the 3.5 m WIYN telescope at KPNO in December 1995 and March 
1996 A 2048 x 2048 S2KB CCD with a plate scale of 0.2'7pixel was used, giving a field of view 
of 6.8'x 6.8' (69 kpc x 69 kpc). Exposures were taken in on-line and off-line filters centered on 
6618 A and 6709 A for the redshifted Hq+[N II] AA 6563 6583 A emission fines, and on 5055 A 
and 5176 A for the redshifted [O III] AA 4959 5007 A emission Unes. The bandwidths of the 
four filters were 74, 71, 32, and 52 A respectively. The average seeing in the images was 0.8" 
except for the B and [O III] images where it was ~ 1.3". We used the IRAF package to obtain 
continuum- free [O III] and Ha-|-[N II] (thereafter referred to as Ha) images. The registration 
accuracy of the optical images is ~ 1". 

The full field of view WIYN R image (Fig. 1, Plate LI) shows the two stellar tails, the 
optical disk, and the three ripples at radii of 25", 45", and 60". Except for these three ripples, 
the optical disk looks relatively undisturbed within a radius of 1' (10 kpc). A magnified view of 
the inner disk region, as seen in the B-V image, is shown in Fig. 2a (Plate L2). It reveals two 
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straight dust lanes which are offset from the nucleus and extend out to the inner ripple. The 
northern dust lane is redder and more prominent than the southern one. The dust lanes bracket 
the circumnuclear region of star formation which is shown in Ha in Fig. 2b (Plate L2). This Ha 
map has a higher resolution (0.8") than previously published maps and resolves the region of 
star formation into (i) a centrally concentrated Ha peak, (ii) a clumpy, arc-like feature about 3" 
(0.5 kpc) north of the peak, extending east-west over 10", and (iii) an Ha bubble with a shell 
feature ~ 6" S of the nucleus. 

In order to unravel the interplay between the molecular gas, the starburst, and the nuclear 
outflow, wc compare this Ha image with the [O III] image, the CO data presented in Paper II, 
and the 5 GHz RC maps from Saikia et al. (1994). This striking comparison is shown in Figures 
3a-3d (Plate L3). The RC map shows a compact central component which is surrounded by 
secondary peaks to the east and west, and two extended bubbles to the north and south. The 
central RC, Ha, and [O III] peaks are coincident within the uncertainties, and lie between the 
two CO lobes (Fig. 3a). The optical emission lines in the inner 3" radius show narrow line 
widths (full width half maximum (FWHM) ^ 300 km s~^) and line ratios typical of HII regions 
(Boer et al. 1992). Thus, the central RC, Ha, and [O III] peaks are likely produced by massive 
star formation. The secondary eastern and western RC peaks lie within the CO bar (Fig. 3a) in 
a region which is associated with a large amount of molecular gas, shows strong Ha emission, 
and has optical line ratios consistent with star formation (Boer et al. 1992). Therefore, the 
secondary RC peaks are also probably due to massive stars. 

However, the northern and southern RC bubbles may have a different origin. They are not 
associated with much molecular gas (Fig. 3a) and are elongated along the CO kinematic minor 
axis which has a position angle of 165 deg (Paper II). The northern RC bubble has a more 
fragmented morphology than the southern one and docs not have optical counterparts. The 
southern RC bubble is associated with the Ha bubble (Fig. 3c) and a bright [O III] bubble (Fig. 
3d), and all three bubbles show a shell feature at 6" (1 kpc) south. At this position, optical 
spectra (Boer ct al. 1992) show that the [O III] lines have a strong extranuclear peak, a large 
FWHM of ~ 580 km s'^, a blueshift of 120 km s'^, and an [O III]/H/3 ratio of ~ 4. The [O III] 
image in Fig. 3d (Plate L3) also shows an extranuclear peak which lies in the shell and has an 
intensity about half that of the nuclear peak. We thus propose that the southern Ha, [O III] , 
and RC bubbles, as well as the northern RC bubble, are part of a starburst-driven outflow 
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emanating from the nucleus. The absence of an optical counterpart to the northern RC bubble 
can be explained in terms of extinction by dust in the optical disk of NGC 2782. The relatively 
undisturbed and circular appearance of the disk (Fig. 1, Plate LI), as well as kinematic evidence 
presented in Paper II, suggests a low inclination for the disk. Furthermore, the near side of the 
disk is probably to the north since the dust lanes are redder and more prominent (Mihalas & 
Binney 1981) in the northern side of the disk (see Fig. 2a, plate L2). We therefore conclude 
that the northern outflow bubble lies behind the disk and is probably obscured from view at 
optical wavelengths by dust in the disk. 

4. The evolutionary phase of the starburst-driven outflow 

The starburst-driven outflow in NGC 2782 has a well-defined collimated bubble morphology 
with a closed shell at its outer edge, as seen in Ha, [O III], and RC. Such an unusual outflow 
morphology has not been previously observed in any other starburst galaxy of comparable 
luminosity (see e.g., Heckman et al. 1993; Lehnert & Heckman 1995, & references therein). In 
most outflows, the Ha and non-thermal RC form complex sets of filaments and loops (e.g., M82, 
NGC 253, NGC 3628, NGC 4666; see Heckman et al. 1993; Seaquist & Odegard 1991; Bland & 
Tully 1988; Fabbiano ct al. 1992; Dahlcm ct al. 1997, & references therein). The soft X-ray and 
RC emission often extends further out than the Ha. The known outflow most similar to NGC 
2782 is that in the LINER galaxy NGC 3079. This galaxy has an Ha bubble of radius 13" (1.1 
kpc), with a shell feature at its outer edge. However, in contrast to NGC 2782, the RC emission 
in NGC 3079 forms two highly fragmented halos which extend twice as far as the Ha bubble 
(Ford et al. 1986; Veilleux et al. 1994). 

What accounts for the remarkable outflow morphology in NGC 2782? Theory predicts that 
in the early stages of a starburst-driven outflow, supernovae and winds from massive stars inject 
energy into the ISM and produce a bubble of very hot gas and thermalised ejecta (e.g., Chevalier 
& Clegg 1985; Tomisaka & Ikeuchi 1988; Heckman, Armus, & Miley 1990). The bubble expands, 
sweeping up ambient gas into a thin, dense, optically-emitting shell. Non-thermal RC can be 
produced when relativistic plasma which is driven along the pressure gradient of the ISM, 
undergoes synchroton losses and inverse Compton scattering (e.g, Heckman et al. 1993; Baum 
et al. 1993; Colbert et al. 1996). Once the bubble is several times the scale height of the 
ambient gas (blowout phase), the swept-up shell fragments through thermal instabilities and 
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allows the hot gas to expand freely. For instance, the starburst galaxy M82 is likely to be in 
this post-blowout, freely-expanding phase (e.g., Heckman et al. 1990), where the dense shell has 
ruptured, allowing the RC-emitting plasma and the hot X-ray-emitting gas to escape beyond 
the optical nebula. On the other hand, NGC 3079 might be in the blowout, partially-ruptured 
bubble phase (Veilleux et al. 1994). We suggest that the southern starburst- driven outflow 
in NGC 2782 is dynamically younger than these outflows; it is probably in a pre-blowout or 
early stages of a blowout phase where thermal instabilities have not yet completely ruptured 
the optical bubble. The electron densities (ne = 600 ± 200 cm^^ ^^^-^ ^^^^ jg 6716 A/[S 
II] A 6730 A line ratio; Boer et al. 1992) in the optical shell of NGC 2782 arc higher than 
densities observed in the outflow region of many galaxies including M82, NGC 253, NGC 3079, 
NGC 3628, NGC 4527, NGC 4536, NGC 4666, NGC 7541, etc (e.g., see Lehnert & Heckman 
1996). This suggests that the shell feature in NGC 2782 contains highly compressed ionised gas. 
The likely enhancement of the magnetic field strength in the compressed material could explain 
the strong RC emission in the shell. 

Our results on the outflow morphology at Ha, [O HI], and RC, together with the short 
outflow dynamical timescale (~ 4 x 10^ years) found from optical spectra by Boer et al. 
(1992), suggest that the outflow is dynamically young. This young dynamical phase could be 
attributed to several factors. The circumnuclear starburst in NGC 2782 may be younger than 
other starbursts which have a comparable luminosity, but a biconical freely-expanding outflow, 
e.g., M82. The dynamical timescale of the outflow in NGC 2782 is less than in M82, where 
we computed the timescale to be ~ 10'' years, using the size and expansion velocity of the 
optical nebula (e.g., Heckman et al. 1990). Furthermore, molecular gas is present around the 
circumnuclear starburst of NGC 2782 (see Fig. 3a, Plate L3) and it can remove energy from the 
outflow and slow down its expansion, for instance, by increasing the rate of radiative cooling 
(Heckman ct al. 1993). In addition, the presence of clouds can reduce the density contrast 
between the dense compressed shell and the hot fluid it encloses, thereby delaying the onset of 
thermal instabilities which rupture the shell. 

Knowledge of the distribution of hot ionised gas in the outflow would help assess its 
dynamical phase. In Fig. 4a, we give the Ha map of a larger region of NGC 2782, and in Fig. 
4b, we compare this map with an archival ROSAT HRI X-ray map which has a resolution of 
~ 12". The X-ray map shows a bright core surrounded by fainter emission extending out to a 
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radius of 20" (3.5 kpc). Unfortunately, the low resolution of the X-ray map does not resolve the 
outflow bubbles, and makes it difficult to separate hot gas driven out by the current outflow 
and that blown out by previous episodes of star formation. Such past episodes are not unlikely. 
The Ha image shows diffuse Ha arcs, loops, and filaments around the central region of star 
formation (Fig. 4a, Plate L4), and this morphology is highly suggestive of gas blown out in 
energetic episodes associated with previous starburst-driven outflows or with the interaction. 
The timcscalcs involved also leave open the possibility of recurrent starbursts; the current 
outflow has a dynamical timcscalc of only ^ 4 x 10® years while the time of closest approach for 
the interaction is 2 x 10^ years (Smith 1994). Although the X-ray map has a low resolution, 
it nonetheless suggests an interesting possibility. The X-ray emission is more elongated to the 
north than to the south, and the northern RC bubble is also more extended and fragmented 
than the southern RC bubble (Fig. 3a, Plate L3); together, these facts suggest that the northern 
part of the outflow may have already broken out and allowed the hot gas to escape. It is unlikely 
that the more extended northern X-ray emission is an instrumental artifact because a similar 
extension is not seen in a point-like source near NGC 2782 in the ROSAT map. Furthermore, 
it is striking that the X-ray emission appears to extend out to the northwestern Ha arc whose 
narrow, bright appearance is suggestive of highly shocked gas (see Fig. 4b, Plate L4). 



5. Cold, warm and hot gas in the outflow 

Does the nuclear outflow in NGC 2782 contain cold molecular gas in addition to hot and 
warm ionised gas? Fig. 3a (Plate L3) shows two CO spurs (labelled 01 and 02) which are 
elongated along the kinematic minor axis, and appear to originate from the center of the CO bar 
where the starburst intensity peaks. The CO isovelocity contours show clear kinks at the base 
of these spurs, indicating deviations from circular motions (see Paper H). The spatial velocity 
plot (Fig. 5 of Paper IT) shows non-circular motions of ~ 30 km s~^ in the two CO spurs. With 
the near side of the disk being north, these velocities are consistent with vertical outflow of gas 
out of the galaxy plane and/or radial inflows of gas in the galaxy plane. The fact that the spurs 
are elongated along the minor axis and lie inside the RC outflow bubbles (Fig. 3a, Plate L3) 
strongly suggests that the spurs contain outflowing gas. 

In order to determine the power needed to drive the outflow, we estimate the energy 
associated with the different phases of the ISM. With a Milky Way value for the CO-H2 
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conversion factor, a radial vertical outflow velocity of (30 km s~^/Sin i ), an inclination i of 30 
deg (see paper II), the CO spurs have ~ 2 x lO'' Mq of H2 and a kinetic energy of ~ 8 x 10^^ 
erg. The warm ionised gas in the optical nebula has a mass of ~ 1.2 x IO^Mq, a kinetic energy 
of ~ 7 X 10^^ erg (Boer et al. 1992), and a thermal pressure of ~ 3 x 10~^ dynes cm~^ for 
600 cm~^ (see § 4). Using the X-ray luminosity of 5 x 10^*^ erg s^^ (from fluxes in Kinney et al. 
1984 for a spectral index of 0) and following Nulsen, Stewart and Fabian (1984), we find that 
the hot gas has a total mass of '--^ 6 x lO'^ f-"'' Mq, a thermal energy of ^ 2 x 10"''^ f*^"^ erg, and a 
pressure of ^ 3 x 10^^" £-0.5 jyj^Qg cm^^. We assumed a temperature of 3 x 10^ K, a radius of 
15" (2.5 kpc) and a volume filling factor f for the hot gas. The above estimates suggest that cold 
molecular gas and hot ionised gas dominate the mass and energy in the outflow of NGC 2782. 
Although the outflows in M82 and NGC 2782 appear to be in different dynamical phases, they 
contain similar masses of cold, warm, and hot gas. In M82, these masses are estimated to be ~ 
a few xlO'^ Mq (Stark & Carlson 1984), 2 x 10^ Mq (Heckman et al. 1990) and 1 x lO'^ Mq 
(Schaaf et al. 1989) respectively. 

The optical linewidths and line ratios in the central region of NGC 2782 do not indicate 
the presence of an AGN (see §3). Furthermore, NGC 2782 has a FIR luminosity (~ lO''^ erg 
s~^) and warm FIR colors (Seo/im/Sioo/im ~ 0.6) similar to those of galaxies with the best 
evidence for supernova-driven winds (e.g., Lehnert & Heckman 1996). (It is however noteworthy 
that the ratios of FIR to blue optical luminosity of both NGC 2782 and NGC 3079 are below 
1 rather than > 2 as is the case in most of these galaxies.) We investigate if star formation 
alone can power the outflow. We assume that the total energy of ~ 2 x 10^^ f*^'^ erg in the 
outflow is a fraction x of the total energy generated by supernovae, where x depends on how 
efficiently the kinetic energy of supernovae is thermalised. With a mean energy of 2 x 10^-'^ 
erg per supernova (Mac Low & McCray 1988) and an outflow timescale of 4 x 10^ yrs (Boer 
et al. 1992), the average supernova rate needed to power the outflow is « (1G~^ f "''/a;) yr"-"^. 
Using the FIR luminosity of 2 x 10^" Lq as an estimate for the bolometric luminosity, and 
applying the starburst model of Rieke et al. (1980) gives a supernovae rate of ~ 0.1 yr^^. For a 
volume filling factor f of 0.1, the supernovae rate from the starburst can power the outfiow if x 
> 5%. Stellar winds can also be an important energy source in establishing a wind (Leitherer & 
Heckman 1995). 

NGC 2782 provides an unparalleled opportunity for setting observational constraints on the 
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poorly-studied, early evolutionary phase of a starburst-driven outflow. Further observations 
of the dynamically young outflow in NGC 2782 can extend the study presented in this paper. 
In particular, future high resolution X-ray and near-infrared observations will help to further 
constrain the dynamical phases of both the southern and northern outflows (§ 4), better 
estimate the poorly-known quantities f and x, and study how the properties of the circumnuclear 
starburst influence the outflow it is driving. 
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Fig. 1 (Plate LI)— The WIYN R image of NGC 2782 with a 6.5' (66 kpc) field of view 
reveals the optical disk bracketed by two stellar tails. Except for three ripples at radii of 25", 
45", and 60", the optical disk is relatively undisturbed within a radius of ~ 1' (10 kpc). 

Fig. 2 (Plate L2) — (a) The WIYN B-V image with a 1.6' (16.3 kpc) field of view shows 
two dust lanes which are offset from the nucleus and extend out to the inner ripple at r ~ 25" 
(4.2 kpc). Lighter regions in the image have redder B-V colors and darker regions have bluer 
B-V colors. The northern dust lane is redder and more prominent than the southern one. The 
dust lanes bracket the inner few kpc region which harbors one of the most luminous starbursts 
among nearby (D< 40 Mpc) spirals, (b) The 0.8" resolution WIYN Ha image with a central 20" 
(3.4 kpc) field of view, shows a central Ha peak, an arc of star formation extending east-west, 
and a southern Ha outfiow which has a well-defined bubble morphology with a shell feature at 
its edge. 

Fig. 3 (Plate L3) — The starburst driven outflow: (a) 5 GHz RC (contours, Saikia ct al. 
1994) on the 2.1" x 1.5" CO map (greyscale, from Paper II) for the central 20" (3.4 kpc). It is 
likely that the main RC peak, as well as the secondary eastern and western peaks, are due to 
massive star formation. In contrast, the northern and southern RC bubbles which are elongated 
along the kinematic minor axis, are associated with little molecular gas and are probably part of 
the starburst-driven outflow. The two CO spurs, labelled 01 and 02, are also elongated along 
the kinematic minor axis and lie inside the RC outflow bubbles. The contour levels are 5, 10, 
30, 50, 70, 90, 100 % of the peak flux (3.97 mJy per beam), (b) The 0.8" resolution WIYN 
Ha image shows a region of star formation whose activity peaks between the two CO lobes, 
and a remarkable southern Ha outflow bubble with an outer shell feature at ~ 6" south. The 
cross marks the 5 GHz RC peak. The contour levels are 65, 75, 85, 90, 95, 100 % of the peak 
flux, (c) Ha (contours) on the the 5 GHz RC (greyscale) . The Ha and RC peak are coincident 
within the uncertainties, and lie between the two CO lobes. The Ha shell at ~ 6" (1 kpc) south 
coincides with the shell in the southern RC bubble, (d) The [O III] image (with 1.3" seeing) has 
a similar bubble morphology as the outflowing Ha, and shows a strong extramiclear peak at ~ 6 
" south. This extranuclear peak occurs in the region where long-slit spectra shows a maximum 
in the intensity and FWHM of blueshifted [O III] lines. 

Fig. 4 (Plate L4) — (a) The Ha image with a 76" (13 kpc) field of view. This figure shows a 
string of HII regions which coincides with the first optical ripple at r ~ 25" (4.2 kpc), a bright 
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northwestern Ha arc, and diffuse Ha loops and filaments which surround the central starburst 
region, (b) The ROSAT X-ray map (contour) of resolution ~ 12" on the Ha image (greyscale). 
There is bright emission in the central 8" (1.4 kpc) radius where the optical and RC outflow 
bubbles lie, and fainter emission which extends further out. The faint X-ray emission is more 
extended northwest than southeast, and it stretches out to the prominent northwestern Ha arc. 
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